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Abstract The ANTARES code has been designed for simulation of astrophysical 
flows in a variety of situations, in particular in the context of stellar physics. Here, 
we describe extensions as necessary to model the interaction of pulsation and con- 
vection in classical pulsating stars. These extensions encomprise the introduction of 
a spherical grid, moveable in the radial direction, specific forms of grid-refinement 
and considerations regarding radiative transfer We then present the basic parame- 
ters of the cepheid we study more closely. For that star we provide a short discussion 
of patterns of the H + He I and the He II convection zones and the interaction with 
pulsation seen in the pdV work or atmospheric structures. 



1 Introduction 

Time-dependent, non-linear models of the classical pulsating variables (cepheids, 
RR Lyr) are around since decades. Nearly all those models have treated the issue 
in one spatial dimension (plus time). As is well known, for the stars near the red 
edge the inclusion of convection is mandatory. Given the basic ID computational 
setting, this had to be done using simplyfied models for convection. For a detailed 
discussion of shortcomings of this approach and their consequences see |[T|,||2|- 

We therefore have adapted the ANTARES code, 15], in order to comply with the 
needs for modelling such stars in multidimensions (2D presently). In section |2] we 
address a few technical issues, in section[3]we discuss a specific model. 
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2 Technical Issues 
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The simulations were performed on a stretched, polar and moving grid. We discuss 
the radiative transfer equation (RTE) and restrictions on time stepping in turn. 

The Radiative Transfer Equation. Near the surface, nontrivial energy exchange 
between gas and radiation is included via the RTE with grey opacities, while dif- 
fusion approximation is used for the deeper layers. The RTE in ID is solved along 
single rays via the short characteristics method of Kunasz and Auer (4)1. In 2D either 
12 or 24 ray directions are chosen according to the angular quadrature formulae of 
type A4 or A6 of Carlson |3], and the directions in each quadrant are arranged in a 
triangular pattern. For each ray the points of entrance and exit plus the correspond- 
ing distances are determined. Since the grid moves this has to be redone every step. 
The RTE 

dl 

^Tr-S-I (1) 

is then solved along each ray. This procedure is repeated recursively since after the 
first step one gets just the intensity on a single new point. 

Time Stepping. The time step for our (explicit) method is never limited by the 
classical CFL condition but is due to radiation transport. Using the diffusion ap- 
proximation also near the surface would enforce prohibitively small time steps, 

At,^ffu.,.e - mm [j^^ (- j j ^ mm j . (2) 

Fortunately, the radiative transfer equation allows the time step to stay in contact to 
the time scale for relaxing a temperature perturbation on the scale of the grid size 
by radiation (Spiegel's relation, ||6l), 

^^™^^'"^"(T6^('-f"f ) ')' 
with k = inin(4r 4 ■■) ' ^hich in practice leads to larger time steps. 



3 Results 

3.1 Model settings 

The parameters of our cepheid are Tg// = 5125 K, \ogg — 1.97, L = 913 L,;;, and 
M = 5 Mq. The period is about 4 days, the radius is 26 Gm, of which the outer 
11.3 Gm were modelled. The opening angle used here is either 1° or 10°. For the 
10° model the computational area has 510 radial grid points and 800 lateral points. 
For the 1° model these number are 800 and 300, respectively. Radially the grid 
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is stretched from cell to cell by a factor of 1.011 and 1.07, respectively. The radial 
mesh sizes of the 10° model vary from 0.47 Mm at the top to 124 Mm at the bottom. 
The runs with these opening angles were started from one and the same ID model 
which had been relaxed for more than 100 periods. - The 10° model is wide enough 
to harbour convection cells of the Hell convection zone. Its resolution is insufficient 
to reasonably represent the H + Hel convection zones. This purpose is achieved by 
the 1° model. - We discuss now the Hell and the H + Hel convection zone in turn. 



3.2 The Hell convection zone 

After some time granted to the Hell c.z. for proper development, ten useful periods 
for evaluation are at our disposal. Visualization shows a clear dependence of con- 
vection on phase. Near maximum compression the lower part of descending plumes 
tends to detach from the body of the c.z., and upon expansion the body reorganizes, 
gains strength and new plumes appear. 

Figure [U shows phase dependence for the pdV work of the Hell zone, averaged 
over the ten useful periods. pdV refers to the work done by convection only, without 
pulsational contribution. Besides the varying strength of the work we see also some 
contributions due to plumes plunging into the stable region below. 

These plumes also generate noticeable gravity waves, discernible down to our 
lower boundary. Laterally, their extent coincides with the width of our computa- 
tional domain. An even wider sector may be appropriate. They have not yet reached 
their statistically steady state, keeping growing. 

With such models, it is also possible to test and calibrate among others convection 
models routinely used in modelling of radially pulsating stars. 




Fig. 1 pdV work due to convection, averaged over ten periods as a function of pulsational phase. 
Below the convection zone proper, overshooting plumes also contribute to some extent. 
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3.3 The H+ He I convection zone 

With the 1° model the H + Hel is reasonably resolved and, as a result, much more 
vigorous than in the 10° case. Velocities (pulsation subtracted) easily get supersonic. 
Note the network of shocks above the hydrogen ionization front in figure |2] 




Fig. 2 Upper convection zone upon contraction. Top is to the left. Colours: density of gas. Isolines 
mark the indicated ionization degrees of hydrogen and helium. 

Frequently, shocks move upwards. In one case, this leads to the strong high- 
density region quite at the top. Here, our closed boundary conditions failed. We 
are working on implementing open boundary conditions at the top. - Mass loss 
in cepheids is a topic of present discussion. Results as the present one hint at the 
intriguing possibility that convectively induced shocks may be a cause. - Generally, 
it should be noted that such calculations are very expensive even in 2D and by no 
means all cepheids are within computational reach. 
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